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Abstract. Discrete-time regulatory networks are dynamical systems on directed graphs, with a 
structure inspired on natural systems of interacting units. There is a natural notion of determination 
amongst vertices, which we use to classify the nodes of the network, and to determine what we call 
sets of dominant vertices. In this paper we prove that in the asymptotic regime, the projection of 
the dynamics on a dominant set allows us to determine the state of the whole system at all times. 
We provide an algorithm to find sets of dominant vertices, and we test its accuracy on three families 
of theoretical examples. Then, by using the same algorithm, we study the relation between the 
structure of the underlying network and the corresponding dominant set of vertices. We also present 
a result concerning the inheritability of the dominance between strongly connected networks. 
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1. Introduction 

In this work we consider a class of discrete-time dynamical systems on networks, inspired on natural 
systems of interacting units as for instance the gene regulatory networks. This class of discrete-time 
models was first introduced in [14j . and in [2] their low dimensional representatives were studied in 
full detail. These systems can be thought as a discrete-time alternative to systems of piccewise affinc 
differential equations extensively studied in [4j[6l[T0], and to finite state models, better known as 
logical networks, which have been used in 0171 [13] for instance. In each of these modeling strategics, 
the interacting units have a regular behavior when taken separately, but are capable to generate 
global complex dynamics when arranged in a complex interaction architecture. In the discrete-time 
dynamical systems considered here, all interacting units evolve synchronously at discrete time steps, 
and the level of activity of each one of them varies according a rule involving its neighboring units in 
the network. 

The interaction architecture may allows to reconstruct the dynamics of the whole systems from the 
observation of a small number of units. The record of the activity of those units during the evolution 
contains all the information needed to determine the state of the system at all times. The aim of the 
present work is to give a topological characterization of these determining units, in the framework of 
the discrete-time dynamical systems on networks mentioned above. This characterization is achieved 
via a classification of the vertices of the network, according to a determination relation we define 
below. The main result states that the vertices in the leading class capture the whole dynamics of the 
network, i. c., if two orbits are such that their restriction to the leading vertices approach to each other, 
then they asymptotically coincide at each vertex of the network. On one hand, the vertices in the 
leading class, or dominant vertices as we call them, are sites to record data in order to reconstruct the 
state of the whole system at each time. On the other hand, those vertices can be used as controlling 
sites. If one forces the configuration of the system to follow a particular trajectory at those sites, then 
this forcing will spread through the whole network. 

l 
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The paper is organized as follows. After some preliminary definitions and notations given in the 
next section, we present in Section [3] the main results, as well as an algorithm determining a set of 
dominant vertices. In order to get an insight in the relation between the underlying network and the 
structure of the dominant set of vertices, in Scction|4]we study several families of theoretical examples. 
Then, in Section [5j we perform a numerical study of the dominant set of vertices in Erdos-Rcnyi and 
Barabasi- Albert random networks. We finish with some concluding remarks and general comments. 

2. Preliminaries 

A discrete-time regulatory network is a discrete-time dynamical dynamical system on a network. We 
start with a network Q = (V,A), where V are the vertices and represent interacting units, and the 
arrows A C V x V denote interaction between them. For each v G V, let I(v) := {u G V : (it, v) G A} 
be the set of vertices influencing v, and 0(v) :={m£F: (v, u) G A} the set of vertices influenced by 
v. To each interaction (u, v) G A, we associate a partition V uv of [0, 1], by a finite number of intervals. 
Each atom in V uv corresponds to an activation mode for the interaction (u, v) G A, so that each atom 
in Q v := riiY^) 'Puv * s associate to a total activation state on the vertex v G V. For each unit u G V, 
the activation partitions {V uv : u G I(v)} define a refinement V u '■= \l ve y'Puvi composed by the 
intersections Hugyluv with I uv G V uv . Finally, these refinements give place to the ff=V dimensional 
partition V :— Y[ v V v . We will denote I(x) to the atom of V containing x. 

The collective activation mode of the network remains constant inside each atom of the partition 
V defined above, and the evolution of the system is given by the iteration of the piecewise affinc 
transformation F : [0, l] v — * [0, such that 

(1) x'+i := F(x') = oxt + (1 - o)£>(x*), 

where a £ [0,1) is the contraction rate related to speed of degradation of the activity of the units 
in absence of interaction, and D : [0, 1} V — > [0, 1} V is a piecewise constant function representing the 
regulatory interactions. The piecewise constant term D : [0, l] v — > [0, 1} V depends on the interaction 
architecture, and for each v € V is such that 

(2) £>(x) t , = D v (x u : u G I{v)) = D v (x 7(t ,)) , 

where D v : [0, l] 7 ( u ) — > [0, 1], is constant in the atoms of Q v 0. Finally, the discrete-time regulatory 
network defined by the interaction function D and the contraction rate a, is the discrete-time dynam- 
ical system ([0, F), with F as defined in (Q]). The interaction architecture, codified in the digraph 
(V, A), is implicit in the definition of the interaction function. 

As usual, the orbit corresponding to the initial condition x° G [0, 1] , is the sequence {x'}t 6 N in 
[0, l] v such that x t+1 = F(x') for each t G N. We will also consider controlled trajectories, which we 
define as follows. For U C V fixed, a U -controlled trajectory is any sequence {x'} 46N in [0, l] v such 
that x^ijy = F(x t+1 )y\u. Here, the projection {x£;} t6N , which we also call control term, is taken 
arbitrarily. 

The discontinuity set A-p of the transformation F is composed by the borders of all the rectangles 
in V . It is a union of (#V — l)-dimensional hyperplanes perpendicular to the principal axes of the 
hypercube [0,1]^. We will say that the orbit or controlled trajectory {x*}^ is uniformly separated 

4n previous works [2] 1111 it has been considered less general interaction functions defined as follows. For each 
interaction (u, v) G A there are associated a sign a uv 6 { — 1,1}, and a threshold T uv G (0,1). Only two activation 
modes are considered: activation if a uv = 1, and an inhibition in the case a uv = —1. The interaction function 
D : [0, l] v -> [0, l] v , is then given by D(x)„ = (#7(d))- 1 Y,uev H(a uv (x u - T uv )). with H : R -> K the Heaviside 
function. 
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from the discontinuity set, if there exists e > such that initio (inf zS A P |x* — z|) > e. Following the 
arguments developed in [3], it can be shown that uniform separation from the discontinuity set is a 
typical situation. For small enough contraction rate, it is should be generic (in the topological sense), 
with respect to the choice of the partitions. In general, by using an appropriate distribution in the 
set of partitions, the uniform separation could be proved to hold with probability one. 



3. Dominance 



For vertex sets U, U' C V, we say that U determines U' . if U 2 Uuet/' i- e -i a ^ the arrows with 
head in U' have their tails in U. We will denoted by dU the maximal set determined by U, which is 
given by 

(3) dU := {v e V : I(v) + and I(v) C U}. 

Whenever U determines U', D(x) v at the coordinates v S V can be computed from the values 
x M at coordinates u S U. Indeed, if xjj and yu belong to the same rectangle in Y[ u eu'^' u ^ then 
D(x)w = D(x)w Since Dw = Yl ue u' ® u ( x -f(«))' an< ^ ^ — UveU' I( v )> tnen we wm use *h e notation 
D(x)u r = D(xu) in this case. 

A vertex set U C V will be called dominant if there exists a nested sequences of vertex sets, U = Ui C 
U 2 C • ■ ■ Ud C f/d+i = V, such that 

dUi = U i+1 \Ui, i = l,...,d. 

Clearly U uniquely determines the sequence U2 ■ ■ ■ Ud+i- We will refer to its length d, as the depth 
of the dominant set U . The classical definition of dominance for directed graphs (as given in [8] for 
instance), stipulates that U is dominant if V = U U dU. Here we slightly extend this notion to allow 
determinations in more that one step. 

According to our definition, a dominant set of vertices U C V governs the dynamics of the whole 
network through a chain U = U% C U2 ■ ■ ■ Ud+i — V of set, so that for each i = 1, . . . , d, Ui determines 
Vi := Ui+i \ Ui. Since IJ-uey I( v ) — then D(x) v at the coordinates v G Vi can be computed from 
the values x M at coordinates u £ U. This allows us to establish the following chain of equations, 

(4) F{ XVl ) = ax Vl + (1 - a)D(xu) 

F(xy 2 ) = axy 2 + (1 - a)D(xu 2 ) 



F(x Vd ) = ax Vd + (1 - a)D(x Ud ) 
F(xcj) = axu + (1 - a)D(x Ud+1 ), 

for all x e [0,l] v . 

The equations in ^ are the key ingredient behind the proof of the following theorem. Other im- 
portant clement is the concept of uniform separation which we remind now. The orbit or controlled 
trajectory {x'j^g is uniformly separated from the discontinuity set, if there exists e > such that 
inf t ^o (inf z( =Ap |x* — z|) > e. We have the following. 

Theorem 3.1. Let U C V be a dominant set of vertices, and {x*} te N, {y*}teN two orbits or U - 
controlled trajectories, both uniformly separated from the discontinuity set. If limt_»oc \x\j — y\j\ = 0, 
then lirm^oo |x* — y'| =0. 
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Proof. Let U = Ui C U2 C Ud+i = V be the nested sequence associated to U, and for each 
i = l,2,...,p, let Vi = Ui+i \ Ui. Let e > be such that inf^o (inf ze A P |x* — z|) > e and 
inf t ^o (infzeA-p |y* — z|) > e, and take an arbitrary <5 < e/ (2^/d + 1). 

Starting from time t — to, we iterate the equations in (|4|) in order to obtain, for each tq € N, 



(5) 



,-r+to 



.T + T0+t0 

c v 2 



£(x 4 +<°) 



T-l 

T-l 

^x^ +t « + (1 - a) J2 « T " 1_t £(x£ T ° +t °) 



T-l 



r+(d-l)T +t 



(d-l)TO+t 



(l-ajxy- 1 -*^* 



t+(d-l)T +t 01 



t=0 
,i 



and all r ^ 0. A similar set of equations is satisfied by y*. Since lim^,^ |x^ — y^| = 0, and both 
are uniformly separated from the discontinuity set, then there exists t <G N such that for each t ^ t , 
and y\j belong to the same atom in Y\ ue uV u ^ and such that |x^ — y\j\ < S. Hence we have 
-D(x^) = D(y t u ) for t ^ to, and taking into account the first equation in ([5]) we obtain 



,-r+to 



v T + t I 

JU 2 



r T+t 



r T+t 



y T u +t0 ) 



? (x& - y&) + (1 - a) a^ 1 "' (£>(4 + '°) - ^(y^ 40 )) 



t=o 



y^°)©« T («&-y&)| 



,T+t 



r + io 

y<7 I 



1 | x Vi yyj 



By taking r such that a T ° < S, we ensure that |x^ 2 — y\j 2 1 < < e/2 for each i > to + to- 

Therefore and y^ belong to the same atom in Y\ ue u 2 T , u for all t ^ tq + to, which implies 
£>(xj} 2 ) = D(yjj 2 ) for all r ^ tq + to- By using this, and the second equation in ([5]), we get 



,T+T + t 



JU3 



,T + T +t 

c c/ 2 



„ rT+To+to} 
JU2 ) 



Ta+ta 



,T+TQ + t 



YU: 



T+T +t \ 2 1 „2t |„ T "+*0 _ i? T 0+*0| 2 



Mb 



yy 2 



Therefore, |x^ 3 — y^ | < y/35 < e/2 for each t ^ 2 To + to, which ensures that Z^x^) = D(y\j^) for 
each t ^ 2 To + to- Iterate this argument, we finally obtain 



r+(d-l)T +t 



r + {d-l)T +t 



T+{d-l)T„+t _ 
T+(d-l)T + t 



+ (d-l)T +t 



^+ (d _ l)T0+t0 ) 



fx^" 1 



Ko + to (d-l)T +t 

yv d 



T+(d-l)T + t 

l u d 



T + {d-l)T +t 

y u d 



(d-\)T +t 



(d-l)r +t 



In this way we have prove that |x* — y'| ^ y/d + 1 5, for all t c?to + to, and since <5 > is arbitrary, 
the theorem follows. □ 



According to this theorem, the dominant set of vertices U C V captures the dynamics of the whole 
network Q := (V, A) through the chain U = Ui C U2 ■ ■ ■ Ud+i = V, i. e., in the asymptotic regime, the 
projection of the dynamics on the dominant set allows us to determine the state of the whole system. 
This is achieved with a delay which depends on the depth of the dominant set, the contraction rate, 
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and on the particular trajectory. If state of the dominant set is forced to follow a prescribed trajectory, 
the asymptotic dynamics of the whole system will depend only on this forcing. In this case we say 
that the dynamics of the network have been controlled. 



Dominance is inherited from a completely connected networks to their completely connected subnet- 
works. To be more precise in this respect, let us state the following. 



Proposition 3.1 (Inhering Dominance). Let Q = (V,A) be a strongly connected network, and Q := 
(V, A) =4 Q , a subnetwork strongly connected as well. Let U £V be a dominant set of vertices for Q. 
Then, there exists a set of vertices U £U, which is dominant for Q . 



Proof. For v £ V, let Ig(v) := {u £ V : (u,v) £ ,4} and Ig(v) := {u £ V : (u,v) £ A}. Then, for 
V C V, let d Q V := {v £ V : Ig(v) + and Ig(v) C V}, and similarly for d g V . 

Since A £ A, and Q is strongly connected, then, for each V C V, we have dgV 3 dgV . Indeed, 
v £ dgV if and only if Lg(v) ^ and Lg(v) C V. Since Ig(v) C Ig(v) and Q is strongly connected, 
then Ig(v) ^ and Ig{v) C V as well. Therefore v £ dgV . 



Let U = Ui C Us ■ ■ ■ C Ud+i = V be the nested sequence associated to U . For each i = 1, 2, . . . , d, we 
have Ui+i = UiUdgUi. Hence, by taking U = U\ := U, and for each i = 1, 2, . . . ,p, Ui + i := UiUdgUi, 
we ensure that Ui 3 U for each i = 1, 2, . . . ,p. In this way we associate to U = U a nested sequence 
U = U\ C U 2 ■ ■ ■ C Ug +1 = V, with d := min{l ^ i ^ d : U i+ i = V}, which satisfies the requirements 
for U to be a dominant set. □ 



3.1. Determining a Dominant Set. 



The following is an algorithm we have implemented for the determination of a dominant set of vertices. 
Here we use the definition in Equation |3]). 
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Algorithm Dominant Vertices. 




Procedure Initial-Solution. 




Input g := (V, A) . Output : U and U' = V \ U . 




Set U := {u 6 V : Id(w) =0}, and U' := 0. 




For i = l to #1/, W := [/ U t/' , [/' = <9W, End. 




If W = V , then finish, else 




Set V' := V\ W. 




While W / V , 




Set U ■— U U {u} , with u 6 V such that 




Id(tt) minimal, given that Od(u) is maximal. 


For i = l to #V, W — UUU', U' = dW, 


End. 


Set V :=V\W. 




End. 




End. 




Procedure Optimization. 




Input : U and U' = V \ U . Output : U and U' = V \ U . 




Fix an order in U , U := {ui, . . . , • 




For i = 1,2, ...,#{/. 




Set t7:=C/\{^}, C7' :=9C7, and V := [/'u{ui}. 




For i = 1 to #V, W :=UUU', U' = dW , End. 




If fCW, then (7 := [7 and U' = U' , End. 




End. 




Procedure Nested. 




Input : U and t/' = V \ U. Output U = Ui <ZU 2 C U d+1 = V . 




Set Ui :=U , W := t/i , and d := 1 . 




While W^V, U d+1 :=U d UdU d , W:=U d+1 , d = d + 


1, End. 



Notice that domination is a property completely determined by the structure of the underlying network 
Q := (V,A), and does not depend on the particular choice of activation partitions V v , or interaction 
function D. 

4. Theoretical Examples 

In order to understand the relation between the structure of the underlying network and the set of 
dominant vertices, we present below some examples where minimal dominant sets can be explicitly 
determined. For those examples we also know, in a explicit way, the output of the algorithm Dominant 
Vertices. We can therefore test the accuracy of the algorithm in this controlled situation. At the end 
of the section we present a example to illustrate Proposition 13 - 1 1 concerning how the dominance is 
inherited to completely connected subnetworks. 

4.1. Dominant Set of Vertices. 
4.1.1. Regular Trees. 

A b-regular tree with L levels can be constructed as follows. From a finite set B on b symbols, we define 
the vertex set V := L)f~QB e , where B° = {e} with e denoting the empty word. The set B l defines the 
^-th level of the tree, and the arrows are given by A := {(ui, ui) G V : ui G {oj} x B or w € {lj} x B}. 
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In this way, vertices in the £-th level form arrows only with vertices in the adjacent levels B izizl . 
Because of this, the union of odd levels determine the even ones, and vice versa, i. e. 

U # = a(\J iA 

t odd \t even / 

U B < = »(u 

t even \l odd / 

Hence, any of V even : = C vcnB e or V^dd := odd-S^ is a dominant set of vertices. 

Figure Q] shows a 2-regular tree with 3 levels. In this case we can select as dominant set the first level 
in the tree, which has 2 vertices only. 




Figure 1. The 2-regular tree with 3 levels. The dominant vertices are colored in grey. 



If we apply Procedure Initial-Solution to the 2-regular tree with 3 levels, we obtain U = {1, 2}, 
and U' = {£,11,12,21,22}. Procedure Optimization returns the same output. Finally, Procedure 
Nested gives d = 1 and Ui = {1,2} C U2 = V. In general, for a &-regular tree with L level, the 
algorithm Dominant-Vertices returns d — 1 and XJ\ := V dd C U% = V. It is not hard to show that 
V ot id is the smallest dominant set for trees with an odd number of levels. 



4.1.2. Dictatorial Networks. 

In this paragraph wc present a collection of completely connected networks with a single dominant 
vertex. For n 2, consider the networks with vertex set V := {1, 2, . . . , n}, and arrows defined by 
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1 1 1\ 









1 1 1 ... 10 
\l 1 1 ... 11 0/ 

The resulting digraph, which we call dictatorial network on n vertices, was a total of (n— l)(n + 2)/2 
arrows. Figure [2] shows the dictatorial network on 5 vertices, which allows us to understand why 
vertex number 1 is dominant in all dictatorial networks. Here, according to Equation (|3|), vertex 1 
determines vertex 5, then {1,5} determines vertex number 4. The set {1,4,5} determines vertex 3, 
and finally {1, 3, 4, 5} determines vertex number 2. 




Figure 2. Dictatorial network on 5 vertices. 

If we apply Procedure Initial-Solution to the dictatorial networks in Figure [21 we obtain U = 
{1,2,3,4,5} and U' = 0. Procedure Optimization returns U = {1} and U' = {2,3,4,5}. Finally, 
Procedure Nested gives d = 4 and U x := {1} C U 2 ■= {1,5} C U 3 := {1,4,5} C U A := {1,3,4,5} C 
Us := {1, 2, 3, 4, 5} = V. In general, since the dictatorial network on n vertices is such that 

d{l} = {n} 

d{l,n} = {n-1} 

d{l,n-l,n} = {n-2} 



the adjacency matrix 



(0 


1 


1 


1 








1 


1 





1 


1 


1 



5{l,3,...,n} = {2}, 
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then algorithm Dominant- Vertices returns d = n — 1 and Ui = {1} C Ui = {1,^} C U3 = {l,n — 
1, n} ■ ■ ■ C U n = V. Of course, the dominant set that is produced by the algorithm is of minimal 
cardinality. 



4.1.3. Democratic Networks. 



In contrast to the previous family of examples, here we present a collection of networks on 2n vertices, 
whose dominant sets have cardinality of at least n. They are bi-directional circuits on an even 
number of vertices, which we define as follows. The vertex set is V := {0, 2, . . . , 2n — 1}, and arrows 



A := {(u, v) G V x V : u - v 
adjacency matrix has the form 



M A 



1( mod 2n) or u 



/o 1 
10 10 
10 1 
10 


\1 



2n — 1( mod 2n)}. The corresponding 



1\ 









1 1 
1 0/ 



We will refer to the network just defined as the democratic network on 2n vertices. 

Figure [3] shows the democratic network on 6 vertices, the set {0,3,4} determines {1,3,5}, and vice 




Figure 3. Democratic network on 6 vertices. 



Applying Procedure Initial-Solution to the democratic networks in Figure [3 we obtain U = 
{0, 1, 2, 3, 4} and U' = 0. Procedure Optimization returns U = {0, 2, 4} and U 1 = {1, 3, 5}. Finally, 
Procedure Nested gives d = 1 and U\ := {0, 2, 4} C Hi := {0, 1, 2, 3, 4, 5} = V. In general, algorithm 
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Dominant-Vertices returns d = 1 and Ui = {0, 2,4, . . . , In - 2} C U 2 = {0, 1, 2, 3, . . . , 2n - 1} = V. 
This is because 

<9{0,2,4,6,...,2n-2} = {1,3,5, ... ,2n - 1} 
9{l,3,5...,2n- 1} = {0,2,4,6,...,2n-2}. 

Notice that the dominant set that is produced by the algorithm is of minimal cardinality. 

In the three examples examined so far, the dominant set U determined by the Algorithm Dominant- 
Vertices is of minimal cardinality, and the complement U' — V \ U is dominant as well, but this is 
not the general case. Furthermore, even when U' = V\ U is dominant, the number of steps needed to 
determine U from U' can be different from the number d given by the algorithm. This is in particular 
the case of the dictatorial networks. 

4.2. Inhering Dominance. 

To illustrate how dominance is inherited to subnetworks, we will consider three supergraphs of the 
dictatorial network in Figure [2l which we present in Figure l4~2l Starting with the network at left, we 
obtained the other two by eliminating an increasing number of arrows, keeping a strongly connected 
network each time. For the graph at left, {1, 2, 3} is the dominant set of vertices given by the Algorithm 
Dominant -Vertices. Applying the same algorithm to the network at the center of the figure, we obtain 
{1,2} as the dominant set. Finally, for the dictatorial network on 5 vertices, vertex number 1 is the 
dominant vertex the algorithm returns. 




Figure 4. Supernctworks of the dictatorial network on 5 vertices. Their respective 
dominant vertices are colored in gray. 



5. Numerical Examples 

In this section we consider two families of random networks, and we determine a dominant set of 
vertices by means of the Algorithm Dominant- Vertices. We will characterize the structural complexity 
of the network through the expected size of the dominant set U and the number of steps d needed to 
determine V from U. 
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5.1. Erdos Renyi networks. 

An Erdos-Renyi digraph is constructed as follows. Fix a set of vertices V and p £ (0, 1]. The arrow 
set A C V x V is such that for each u,v 6 V, P{(u,v) £ A} = p. According to this, a given arrow 
(u, v) belongs to the digraph with probability p, and the inclusion of different vertices are independent 
and identically distributed random variables (for more details on random graphs and their dynamical 
construction see [5] for instance). 

For p = 0.1, 0.2, . . . , 1 and #y = 50,100,200, we have generated 100 Erdos-Renyi digraphs, and 
for each one of those 3000 networks, we have computed a dominant set of vertices U and a depth 
d, by using Algorithm Dominant -Vertices. On top of Figure [5] we show the expected number of 
dominant vertices E(#C7) the Algorithm Dominant- Vertices returns, and in error bars the size of 
the fluctuations s(#U) := */E(#Z7 — E(#Z7)) 2 , both as a function of p, and for the three scenarios 
#V" = 50, 100, 200. In the lower frame of the same figure we plot the expected number of steps E(d) 
given by the same algorithm, with error bars proportional to the corresponding standard deviation 
s(d) := ^E(d-E(d)) 2 , both as a function of the connection probability p. 

In Figure [5] we show a realization of an Erdos-Renyi network on 100 vertices, with probability of 
connection p = 0.1. For this example, the Algorithm Dominant -Vertices returns #Z7 = 58 and 
d= 13. 

For the completely connected network on N vertices, any dominant set has cardinality N — 1 and 
depth d = 1. Our numerical results show how the Erdos-Renyi networks approaches this situation 
as the probability of connection goes to 1. More specifically, as p — > 1, the size of the dominant set 
E(#t/) -> N, while E(rf) -> 1. 

5.2. Barabasi Albert Networks. 

We implemented a dynamical construction of scale-free graphs similar to the one proposed by Barabasi 
and Albert. We start, at n = 0, with a kernel graph Qq, which is the complete simple graph on too 
vertices. Then, for each n ^ 1, we add a new vertex v n +\ to the preceding graph Q n :— (E n ,V n ). 
This new vertex form new edges with randomly chosen vertices in V n . A given vertex v € V n has 
probability proportional to its degree to be selected, and we have 

s \ c j? -a. i v.'r+ f \ #{u€V n : {u,v} € E n } 
{v n+ i,v\ G E n with probability p n {v) := — — . 

The graph at the (n + l)-th step is Q n +u w ith V n+ i := V n U {v n +i}, and an enlarged edge set 
E n+ \. This random iteration continues until a predetermined number of vertices is obtained. In 
our dynamical construction it is possible to add several edge at each iteration, in opposition to the 
traditional scheme where only one edge is added at each time step. We finally assign directions to 
the edges in EN-m a , by choosing one of the two possible directions with probability 1/4, and both 
directions with probability 1/2 (for a rigorous study of scale-free and related random graphs, see [5]). 

For #V = 50, 100, 200, and m = 0.06 x #V, 0.10 x #V, 0.14 x #V, . . . , 0.38 x #V, we have generated 
100 Barabasi- Albert networks, and for each one of those 2700 digraphs, we have computed U and d, 
by using Algorithm Dominant -Vertices. Figure [7] is the analogous to Figure [5] above. It shows in the 
upper frame the behavior of E(JfU) — mo, the expected number of dominant vertices exceeding the 
kernel of size, as a function of the relative size of the kernel mo/^V. In error bars we plot the size of 
the corresponding fluctuations s{#U) := ^(#[7 - E(#C/)) 2 . The fr ame at bottom shows the mean 
depth E(d) given by the same algorithm, and error bars proportional to the corresponding standard 
deviation s(#d) := y/E(d - E(rf)) 2 , also as a function of the relative size of the kernel. In both frames 
we compare the three scenarios #V" = 50, 100 and 200. 
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In Figure [8] we show a realization of a Barabasi-Albert network on 100 vertices, with kernel size 
mo = 6. For this example, the Algorithm Dominant-Vertices returns #U = 37 and d = 3. 

We observed that the size of the dominant set for a Barabasi-Albert network increases as the relative 
size of the kernel grows, in a similar way as in the Erdos-Renyi case. Nevertheless, there is an 
important difference. Since the kernel is build from a complete graph, the size of the dominant set 
is expected to be at least as large as the kernel. We observe that number extra vertices needed to 
complete the dominant set decreases with the size of the kernel. At the same time, the depth of the 
dominant set follows a non-monotonic behavior. 

6. Final comments 

Our definition of dominant vertices depends only on the structure of the underlying network, and for 
orbits (or forced trajectories) uniform separated from the discontinuity set, those vertices capture the 
dynamical state of the whole network. In the general case, for orbits accumulating at the discontinuity 
set, another definition of domination should be considered. Nevertheless, since uniform separation 
from the discontinuity set is typical, our algorithm determines observation (or control) nodes, in all 
but a negligible proportion of the cases. 

One could interpret our main result as the existence of a effective network, smaller than the underlying 
network, which could be constructed from the dominant set of vertices. The dynamics on the effective 
network would be equivalent to the original one. Though such a network reduction can be attempted, 
specially in the case of a — (instantaneous degradation), the resulting dynamical system is not 
necessarily simpler that the original one. 

The algorithm we propose is non-deterministic, and several alternative formulations can be considered. 
For instance, in the procedure Initial-Solution one chooses vertices with minimal input degree, 
amongst those of maximal output degree. Alternatively, one could choose vertices maximizing the 
output degree, amongst those of minimal input degree. We have considered this alternative, and 
we have found that in general it produces larger dominant sets. In spite of this, we expect that 
this algorithm can be improved, but that was not our aim in this work. Instead, we have tested 
the efficiency of the algorithm on examples where the dominant sets of vertices can be explicitly 
determined, and we have found optimal or almost optimal results. 

Finally, the properties of dominant sets of vertices can be used to characterize the complexity of a 
given network. We have done this for the Erdos-Renyi and for the Barabasi-Albert ensembles of 
random networks. This topological characterization has to be contrasted with a dynamical one, and 
this is what we intend to do in a future work. 
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Figure 5. On top, the expected number of dominant vertices E(#{/), with error bars 
corresponding to the fluctuation s(#U) := \/E(#U - E(#L/)) 2 , both function of 
p. On bottom, the mean depth E(d), with error bars proportional to the corresponding 
standard deviation s(d) := y/E(d — E(d)) 2 , both as a function of the connection 
probability p. In both frames we compare the three scenarios #V = 50, 100 and 200, 
and their corresponding curves appeared ordered from bottom to top. 




Figure 6. Realization of an Erdos-Renyi Network on 100 vertices. The dominant 
vertices determined by the Algorithm Dominant- Vertices are marked with a x . 
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Figure 7. In top, the behavior of E(#J7) — mo as a function of the relative size of 
the kernel m /#V. The fluctuations s(#J7) := \/E{#U - E(#[/)) 2 are indicated 
by the error bars. At bottom, the mean depth E(gQ, and error bars proportional 
to the corresponding standard deviation s(#d) := -\/E(<i — E(<i)) 2 , both as functions 
of the relative size of the kernel. In both frames we compare the three scenarios 
#V = 50, 100 and 200, and their corresponding curves appeared ordered from bottom 
to top. 




Figure 8. Realization of a Barabasi- Albert Network on 100 vertices. The dominant 
vertices determined by the Algorithm Dominant- Vertices are marked with a x . 



